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We discuss a phenomenological model that describes results on diffractive pp and pp cross sections 
and event final states up to the Fermilab Tevatron energy of yfs = 1.96 TeV and use it to make 
predictions for Large Hadron Collider (LHC) energies up to yfs = 14 TeV and asymptotically as 
\fs — > oo. The model is anchored in a saturation effect observed in single diffraction dissociation that 
explains quantitatively the factorization breaking observed in soft and hard pp and pp diffractive 
processes and in diffractive photoproduction and low Q 2 deep inelastic scattering. 
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I. INTRODUCTION 

As we entered a new energy frontier at the Large 
Hadron Collider (LHC) with data collected at y/s = 
900 GeV, 2360 GeV, and 7 TeV from Fall 2009 to Spring 
2010, it became painfully clear that the Monte Carlo 
(MC) simulations designed to represent the collective 
knowledge of the field on diffractive cross sections and 
event final states did not meet the challenge presented to 
them in this new higher energy environment. The most 
commonly used event generators, PYTHIA [l[ and pho- 
JET , were found to disagree not only with the data but 
also with each other. The latter clearly meant that the 
two simulations could not both be right. Therefore, an 
update of the MCs was urgently needed. Because of the 
importance of Minimum-Bias (MB) MC simulations in 
estimating trigger rates, backgrounds, and the machine 
luminosity at the LHC, a "diffraction" workshop was 
organized at CERN on 7 May 2010 0] that brought ex- 
perimentalists and theorists together to exchange ideas 
with the goal of producing a reliable MC generator for 
the LHC. This paper is based on a talk I presented at 
that meeting and an expanded version presented at this 
workshop. 

Diffraction dissociation in pp/pp interactions may be 
defined by the signature of one or more "large" and char- 
acteristically not exponentially suppressed [4J rapidity 
gaps (regions of rapidity devoid of particles) @ in the 
final state. The rapidity gap is presumed to be due 
to the exchange of a strongly-interacting color singlet 
quark/gluon combination with the quantum numbers of 
the vacuum, traditionally referred to as "Pomeron" (IP) . 
Diffractive processes are classified as single diffraction 
(SD), double Pomeron exchange (DPE), also referred 
to as central dissociation (CD), and double diffraction 
(DD). In pp SDp (SD p ), the p(p) dissociates while the 
p(p) remains intact escaping the collision with momen- 
tum close to that of the original beam momentum and 
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FIG. 1: Non-diffractive and diffractive pp interactions. 

separated from the p (p) dissociation products by a for- 
ward gap; in DPE both the p and the p escape, result- 
ing in two forward gaps; and in DD a central gap is 
formed while both the p and p dissociate. The above ba- 
sic diffractive processes are listed below, along with two 
additional two 2-gap processes which are combinations of 
SD and DD and are indicated as SDD: 

TABLE I: Diffractive cross sections, 
acronym basic diffractive processes 
SDp pp -> p + gap + [p -> X p ], 
SD P pp -> [p -> X p ] + gap + p, 
DD pp [p ->• Xp] + gap + [p -> X p ] , 
DPE pp — s> p + gap + X c + gap + p, 

2-gap combinations of SD and DD 
SDDp pp- 
SDDp pp ■ 



V + gap + X c + gap + [p -> X p ], 
[p -> Vp]gap + X c + gap + p. 



Here, 



X p and X c represent clusters of particles in 
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rapidity regions not occupied by the gap(s). The 2-gap 
processes are examples of multi-gap diffraction, a term 
coined by this author to represent events with multiple 
diffractive rapidity gaps. A special case of DPE is exclu- 
sive production, where a particle state is centrally pro- 
duced, as for example a dijet system or a Z boson. 

Below, in Sec. [TT1 (strategy) we outline the method we 
follow to implement an algorithm for a MC simulation, in 
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FIG. 2: Average multiplicity dN/drj (vertical axis) vs. r\ (horizontal axis) for a process with four rapidity gaps, Ar)i(i = 1 — 4). 



Sec. |TTT] (cross sections and final states) we present 
excerpts from previous papers on total and differential 
diffractive cross sections and final states, and in Sec. IIVI 
we conclude. 



II. STRATEGY 

A phenomenology that is used to make predictions for 
the LHC and beyond should be based on cross sections 
and final states that incorporate the current knowledge 
in the field molded into a form that can be extrapolated 
to higher energies. The issues to be addressed is how 
to take into account saturation effects that suppress 
cross sections, and what formulas to use for event final 
state multiplicity, pseudorapidity, and transverse energy 
(Et) distributions. In addition, the structure of an 
algorithm for implementing this knowledge into a MC 
simulation should also be addressed. The algorithm 
must be robust against changes in the collision energy, 
so that it may be equally well applied to simulate 
collisions at fixed target energies as well as at the higher 
energies of pp and pp colliders and in astrophysics. In 
this section, we outline a strategy that addresses these 
issues. 

The following input is used for cross sections and final 
states: 

(i) d 2 a I d£dt of the diffractive processes listed in TableU 
from the RENORM model 0] ; 

(ii) <Jt(s) from SUPERBALL MODEL 

(iii) optical theorem — > Imf e i(t = 0) (imaginary part 
of the forward scattering amplitude); 

(iv) dispersion relations — > Re f e i(t = 0), using low 
energy cross sections from global fit [9(; 

(v) final states: use "nesting" to describe gap processes, 
where a nest is defined as a region of A77 where there 
is particle production, in contrast to a gap region 
where there are no particles [l(| [H| • 

Figure [2] shows a schematic rj topology of an event with 
four rapidity gaps and three nests of final state particles. 



The cross section for this configuration is presented in 
Sec. InTAl 

We propose the following algorithm for generating final 
states: 

• start with a pp — > X inelastic collision at y/s; 

• decide whether the collision is ND or diffractive 
based on the expected cross sections; if ND, use 
the ND final state expected at s/s; if diffractive, 
select SDp, SD p , DD, or DPE based on probabili- 
ties scaled to the corresponding cross sections; 

• for each diffractive event, check whether the region 
of 77 where particles are produced, A?/, is large 
enough to accommodate additional diffractive ra- 
pidity gaps: if yes, decide whether or not the event 
will have other gaps within this region, again us- 
ing probabilities scaled to the cross sections, and 
branch off accordingly; 

• continue this process until the region Arj' is too 
small to accommodate another diffractive gap. 

It is important to note that in our definition of a ND 
collision there are no diffractive gaps whatsoever in the 
final state of the event. In this respect, this definition dif- 
fers from those of "inclusive" or "non-single-diffractive" 
definitions of ND events used in the literature. 



III. CROSS SECTIONS AND FINAL STATES 

In this section, we discuss briefly the diffraction disso- 
ciation and total cross sections using information and / or 
excerpts from Refs. @, EH- 



A. Diffractive cross sections 

In Ref . [IH , the following expression is obtained for the 
SD cross section [quoting] : 
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where: 

(i) the factor in square brackets represents the 
cross section due to the wee partons in the rj- 
region of particle production At/; 

(ii) Arj — In s-Arj 1 is the rapidity gap; 

(hi) k is a QCD color factor selecting color- 
singlet gg or qq exchanges to form the rapidity 
gap; 

(iv) P gap (Ar/,t) is a gap probability factor 
representing the elastic scattering between 
the dissociated proton (cluster of dissociation 
particles) and the surviving proton; 

(v) Ng a p(s) is the integral of the gap prob- 
ability distribution over all phase-space in t 
and At;; 

(vi) F p (t) in Pg ap (Ai],t) is the proton form 

factor F p (t) = e bot ... ; and 

(vii) Cg a p is a normalization constant, whose 
value is rendered irrelevant by the renormal- 
ization division by N gap (s). ... 

By a change of variables from Aij to M 2 using 
At/ = InAf 2 and At? = Ins - InM 2 , Eq. {J) 
takes the form: 
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where b = b + 2a' In ™ [b is the slope of the 
diffractive ^-distribution] . Integrating this ex- 
pression over M 2 and t yields the total single 
diffractive cross section, 
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constant = a. 



(1) 



intercept to the slope of the Pomeron trajectory. Good 
agreement with the ratio extracted from measurements 
is obtained, providing support for the renormalization 
approach used in the phenomenology. 

A similar expression may be use for DD, DPE, and 
multigap processes, as discussed in Refs. El- For 
example, the differential cross section for the process dis- 
played in Fif. [5] is derived in Ref. as [quoting] : 



(2) 



(3) 



The remarkable property that the total single 
diffractive cross section becomes constant as 
s —} oo is a direct consequence of the coher- 
ence condition required for the recoil proton 
to escape the interaction intact. This condi- 
tion selects one out of several available wee 
partons to provide a color-shield to the ex- 
change and enable the formation of a diffrac- 
tive rapidity gap. 

Details are presented in Ref. [llj], where this formula- 
tion of the cross section is used to derive the ratio of the 
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where the term in square brackets is the pp 
total cross section at the reduced s-value, de- 
fined through ln(s'/so) = Yli K ( one f° r 
each gap) is the QCD color factor for gap for- 
mation, the gap probability is the amplitude 
squared for elastic scattering between two 
diffractive clusters or between a diffractive 
cluster and a surviving proton with form fac- 
tor F 2 (t), and N gap is the (re)normalization 
factor defined as the gap probability inte- 
grated over all 10 independent variables U, 
■q^ril, and Ar] = J2t=i A Vi- 

The renormalization factor N gap is a func- 
tion of s only. The color factors are c g = 
(N 2 — 1) _1 and c q = 1/N C for gluon and quark 
color-singlet exchange, respectively. Since the 
reduced energy cross section is properly nor- 
malized, the gap probability is (re)normalized 
to unity. The quark to gluon fraction, and 
thereby the Pomeron intercept parameter e 
may be obtained from the inclusive parton 
distribution functions (PDFs) |l3j. Thus, 
normalized differential multigap cross sec- 
tions at t = may be fully derived from in- 
clusive PDFs and QCD color factors without 
any free parameters. 

The exponential dependence of the cross sec- 
tion on Arji leads to a renormalization fac- 
tor ~ s 2e independent of the number of gaps 
in the process. This remarkable property of 
the renormalization model, which was con- 
firmed in two-gap to one-gap cross section ra- 
tios measured by the CDF Collaboration (see 
Ref. [lj]), suggests that multigap diffraction 
can be used as a tool for exploring the QCD 
aspects of diffraction in an environment free 
of rapidity gap suppression effects. The LHC 
with its large rapidity coverage provides the 
ideal arena for such studies. 



(4) 
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• Use m 2 = s a in the Froissart formula multi- 
plied by 1/0.389 to convert it to mb -1 . 

• Note that contributions from Reggeon ex- 
changes at yfs = 1800 GeV are negligible.as 
can be verified from the global fit of Ref. [9j . 

• Obtain the total cross section at the LHC: 

„ / „LHC CDF \ 
LHC = cdf + JL . ( In 2 £ In 2 £ ) 

S \ SF SF ) 

For a numerical evaluation of <j LHC we use 
as input the CDF cross section at y/s = 
1800 GeV, erf DF = 80.03 ± 2.24 mb, the 
Froissart saturation energy y/s F — 22 GeV, 
and the parameter s a . 

...The resulting prediction for the total cross 
section at the LHC at y/s = 14 TeV is: 

ofF&s = ( 80 ± 3) + (29 ± 12) = 109 ± 12 mb. 
For y/s = 7 TeV, the predicted cross section is: 
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FIG. 3: Total pp/pp single-diffraction dissociation cross sec- 
tion data (sum of both p and p dissociation) for £ < 0.05 
compared with predictions based on the standard and the 
renormalized Pomeron flux (from Ref. 0). 



B. The total cross section 

In Ref. Q, an analytic expression is obtained for the 
total cross section using a parton model approach and 
exploiting a saturation effect observed in the SD cross 
section. The abstract of Ref. Q reads [quoting] : 

The single-difitactive and total pp cross sec- 
tions at the LHC are predicted in a phe- 
nomenological approach that obeys all uni- 
tarity constraints. The approach is based 
on the renormalization model of diffraction 
and a saturated Froissart bound for the total 
cross section yielding a t — (tt/s ) • ln 2 (s/sF) 
for s > sf, where the parameters s and sf 
are experimentally determined from the y^s- 
dependence of the single-diffractive cross sec- 
tion. 

The following strategy is used in Ref. |8| [quoting]: 

• Use the Froissart formula as a saturated cross 
section rather than as a bound above Sf- 

a t (s> s F ) = a t (s F ) + % -In 2 £ 

• This formula should be valid above the knee 
in a sd vs. y/s at \fs F = 22 GeV (Fig. [3]) and 
therefore valid at y/s = 1800 GeV. 



o-LHC v = Q8 ± 8 mb r at ^ = 7 TeV j ) 

The result for yfs = 14 TeV is in good agreement with 
o-CMG = H4± 5 mb obtained by the global fit of Ref. 
where the uncertainty was estimated from <5e and the s e 
dependence from which the value of the parameter s was 
obtained. 

IV. CONCLUSIONS 

We briefly discuss a phenomenological model that de- 
scribes available results on diffractive pp and pp cross 
sections and event final states up to the Fermilab Teva- 
tron energy of y/s — 1.96 TeV and refer the reader to 
previous publications for further details. We also outline 
a procedure to be used to implement the predictions of 
the model into a Monte Carlo simulation that is robust 
against changes in the collision energy, so that it may 
be equally well applied to simulate collisions at fixed 
target energies as well as at the higher energies of the 
Tevatron, the LHC, and beyond. The model is anchored 
in a saturation effect observed in single diffraction dis- 
sociation that explains quantitatively the factorization 
breaking observed in soft and hard pp and pp diffractive 
processes and in diffractive photoproduction and low Q 2 
deep inelastic scattering. 
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